Footprints in Sand: The Response of a Granular Material to Local Perturbations 
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, We experimentally determine ensemble-averaged responses of granular packings to point forces, 

and we compare these results to recent models for force propagation in a granular material. We used 
C^) \ 2D granular arrays consisting of photoelastic particles: either disks or pentagons, thus spanning the 

range from ordered to disordered packings. A key finding is that spatial ordering of the particles is 
a key factor in the force response. Ordered packings have a propagative component that does not 
occur in disordered packings. 

PACS numbers: 46.10.+Z, 47.20.-k 

Granular systems have captured much recent interest due to their rich phenomenology, and important applications 
[1]. Even in the absence of strong spatial disorder of the grains, static arrays show inhomogeneous spatial stress 
profiles called stress (or force) chains [2] . Forces are carried primarily by a tenuous network that is a fraction of the 
total number of grains. 

A fundamental unresolved issue concerns how granular materials respond to applied forces, and there are several 
substantially different models. A broad group of conventional continuum models (e.g. elasto-plastic, . . .) posit an 
elastic response for material up to the point of plastic deformation [3] . The stresses in portions of such a system below 
q ' plastic yield have an elastic response and satisfy an elliptic partial differential equation (PDE); those parts that are 
plastically deforming satisfy a hyperbolic PDE. Several fundamentally different models have recently been proposed. 
' The g-model of Coppersmith et al. [4] assumes a regular lattice of grains, and randomness is introduced at the contacts. 
This model successfully predicts the distribution of forces in the large force limit, as verified by several static and 
f***- ' quasistatic experiments and models [4-6]. In the continuum limit, this model reduces to the diffusion equation, since 
the forces effectively propagate by a random walk. Another model (the Oriented Stress Linearity-OSL-model) of 
Bouchaud et al. [7], has a constitutive law, justified through a microscopic model, of the form a zz — fia xz + r/a xx 
(in 2D) in order to close the stress balance conditions duij/dxj — pgi. This leads to wave- like hyperbolic PDEs 
describing the spatial variation of stresses. In later work, these authors considered weak randomness in the lattice 
[8], and proposed a convection-diffusion (C-D) equation. A model proposed by Kenkre et al. [9] combines pure wave 
propagation and diffusion in the telegraph equation. Recently, Luding et al. [10] and Bouchaud et al. [11] have 
proposed models that account for force fluctuations on both the length scale of grains and on the length scale of stress 
y ■ chains. 

, The range of predictions among the models is perhaps best appreciated by noting that the different pictures predict 
qualitatively different PDEs for the variation of stresses within a sample: e.g. for elasto-plastic models an elliptic or 
hyperbolic PDE; for the g-model, a parabolic PDE; and for the OSL model without randomness, a hyperbolic PDE. 
O ■ The impact of equation type extends to the boundary conditions needed to determine a solution: e.g. hyperbolic 

, equations require less boundary information than an elliptic equation. 

Here, we explore these issues through experiments on a 2D granular system consisting of photoelastic (i.e birefringent 
under strain) polymer particles [6] that are either disks or pentagons. By viewing the particles through an arrangement 
of circular polarizers (a polariscope) it is possible to characterize the stress on the particles [6]. In the images shown 
below, bright regions correspond to locally large force/stress. Near a contact, the stresses within a particle are very 
nonuniform. This leads to a series of light and dark bands in the polariscope image intensity, /, with the density of 
bands increasing monotonically with the force at the contact. We exploit this fact to produce a force calibration in 
terms of G 2 = |V/| 2 , since the mean (G 2 ) over the scale of a disk grows with the contact forces. We obtained this 
calibration either by applying known forces F to the boundary of a small number of particles and at the same time 
measuring (G 2 ), or by applying uniform loads to the upper surface of a large rectangular sample. As a confirmation 
that this technique measures the force (pressure), we show the hydrostatic pressure vs. depth for two different samples 
in Fig. 1. 
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We used two different particle types, disks and pentagons, in order to construct packings of varying degrees of 
spatial order/disorder. For monodisperse disks, we obtained a highly ordered packing. With bidisperse distributions, 
we modified the disorder in a controlled way [12] characterized by the parameter A = (a) 2 /(a 2 ). The brackets refer 
to averages over the sample for powers of the disk diameters a. To vary A, we used mixtures of two disk diameters: 
ai w 0.7 cm and ai w 0.9 cm. Both types of disks had a thickness of 0.64 cm. The pentagons had the same thickness 
as the disks, and a side length of 0.7 cm. A sample, typically <~ 60 particles wide and 25 particles high, was placed in 
a nearly vertical plane. To keep the sample from collapsing, it rested against a smooth powder-lubricated glass plate 
that was inclined from the vertical by a small angle, ~ 3°; hence minimizing the friction with the plate. The sample 
rested on a rigid metal base, and was confined at the sides by rigid metal bars. The samples were prepared by gently 
adding particles to the upper surface until the full amount of particles was in place. 

Fig. 2a shows a typical polariscope image in the absence of any applied force, besides gravity. The stress at the 
bottom of the sample is larger than at the top due to gravitational head, and the mean stress was a linearly increasing 
function of depth. 

We measured the system response by placing a known weight carefully on top of one particle in the approximate 
center of the sample, producing a local vertical force (Fig. 2b). We then removed the weight. For the monodisperse 
and pentagonal packings described below, we used only responses where the particles remained undisturbed by this 
process. For polydisperse disk samples, some small rearrangement of the surface grains always occured, however. For 
each realization, we obtained the stress difference between successive images (Fig. 2c), containing only the response 
from the point perturbation, with the linear hydrostatic head effects removed. 

We repeated this process for many different rearrangements of the particles, typically 50 times for each set of 
particles. The images of Fig. 2c, d show that the responses are complicated and differ significantly from realization 
to realization. There are at least two approaches to address the large variability among realizations: 1) obtain 
information on an ensemble of realizations; or 2) perturb the system over a larger number of grains in order to obtain 
a larger-scale averaging. Here, we pursue the first option. 

In Fig. 3, we show using greyscale the average response for monodisperse disks, a mixture of disks, and pentagons, 
for a point load of 50 g. In Fig. 4, we show more quantitative data for the force response along a series of horizontal 
lines at depths, z (measured from the top). These data show clear differences depending on the packing geometry. The 
monodisperse packing shows two peaks emerging from the central source, broadening with z, but remaining clearly 
identifiable. For the bidisperse packing, the peaks are also present, but much less sharply resolved. Data for the 
pentagons show no evidence for wave-like stress propagation. There is some noisy structure for the pentagon response 
at the greatest depths. We believe this is due to the fact that it is relatively difficult to rearrange the interlocked 
pentagonal packing at large depths; hence chains tended to appear there in nearly the same places. In any event, this 
effect only occurs deep in the pile. 

A key question is then which model best describes these results. The answer clearly depends on the amount of 
spatial order in the packing. For the ordered (monodisperse) packing, the C-D model proposed by Claudin et al. [8] 
provides reasonable agreement with the data, as shown by the fit lines of Fig. 4. However, there are clear departures 
between the data and fits to the C-D model, since the central peak persists to depths for which the C-D model predicts 
no such feature, and the two side structures fall to the noise level for h > 8d. The recent models proposed by Luding 
ct al. [10] and by Bouchaud et al. [11] predict the presence of the central feature, and thus offer improved fits (details 
to be discussed elsewhere). For the disordered packings, best characterized by pentagons, the one central peak might 
be described by either diffusion or by elasticity. 

It is possible to distinguish which of these is best from the width, W(z), of the response peak. For a diffusive model, 
W should grow as \fz. For an clastic layer, the response to a point force is a modified Lorentzian with W ct z for 
moderate z. We show in the inset of Fig. 4c data (on log scales) indicating that W oc z over much of the experimental 
range, as in an elastic response. 

The pentagon results, including the linear dependence of W on z, are consistent with recent experimental work 
by Reydellet et al. [13] for 3D disordered packings, with numerical simulations by Moreau [14], and with the recent 
theoretical work of Bouchaud et al. [11], who find a crossover to an elliptic response beyond a length scale controlled 
by the level of disorder. Recent experiments have also been carried out by Rajchenbach [15], who used cuboidal 
blocks of photoelastic material in a similar arrangement to that used here. These experiments showed a parabolic 
outer envelope of the response. Since it is difficult to know the nature of the intergrain contacts and the extent of 
their order/disorder, it is also difficult to place these measurements in the present context. 

To conclude, we have measured the force response of 2D granular packings to local force perturbations. There are 
large variations for any given response, and we consider averages over many realizations. For packings with strong 
spatial order, the mean response has a strong propagative part, and convection-diffusion or more recent models [10,11] 
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give a reasonable description. As the amount of disorder increases, the propagative component diminishes, and the 
response is most like that of an elastic material. Interestingly, the connection between order/disorder and the force 
response may provide a way to characterize the disorder in granular samples, something that may be important for 
understanding experiments where the amount of disorder of a packing was clearly important [16]. These results raise 
additional questions for future work such as the need to investigate the vector character of force propagation (e.g. 
forces applied at arbitrary angles to a surface, forces applied in the interior). Another important issue concerns the 
statistical variability that can be expected in a single realization. 

Acknowledgments We appreciate helpful interactions with P. Claudin, I. Goldhirsch, D. Schaeffer, and J. Socolar. 
The work of SL was supported by the Deutsche Forschungsgemeinschaft (DFG) and the National Science Foundation. 
The work of GR, LV, and EC was supported by PICS-CNRS #563. The work of JG, DH, EL, and RPB was supported 
by the US National Science Foundation under Grant DMR-9802602, and DMS-9803305, and by NASA under Grant 
NAG3-2372. 



[1] For a broad perspectives see Focus Issue on Granular Materials, Chaos 9, 509-696 (1999); Physics of Dry Granular Media, 
H. J. Herrmann, J.-P. Hovi, and S. Luding, eds. NATO ASI Series, Kluwer, 1997; Powders and Grains 97, R. P. Behringer 
and J. T. Jenkins, eds. Balkema, 1997; H. M. Jaeger, S. R. Nagel, and R. P. Behringer, Rev. Mod. Phys. 68, 1259 (1996). 

[2] S. Ouaguenouni and J.-N. Roux, Europhys. Lett, 39 (2), pp. 117-122 (1997) 

[3] D. M. Wood, Soil Behaviour and Critical State Soil Mechanics (Cambridge University, Cambridge, England, 1990); R. M. 

Nedderman, Statics and Kinematics of Granular Materials, Cambridge Univ. Press, Cambridge, 1992. 
[4] C.-h. Liu, S. R. Nagel, D. A. Schecter, S. N. Coppersmith, S. Majumdar, O. Narayan and T. A. Witten, Science 269, 513 

(1995) , Phys. Rev. E 53, 4673 (1996). 

[5] D. K. Mueth et al. Phys. Rev. E 57, 3164 (1998); B. J. Miller, C. O'Hern, and R. P. Behringer, Phys. Rev. Lett. 77, 3110 

(1996) . 

[6] D. Howell, C. Veje, and R. P. Behringer, Phys. Rev. Lett. 82, 5241 (1999); Phys. Rev. E 59, 739 (1999); Chaos 9, 599 
(1999). 

[7] J.-P. Bouchaud, M. E. Cates, and P. Claudin, J. Phys. I, France 5, 639 (1995); J. P. Wittmer, M. E. Cates, P. Claudin, J. 

Phys. R., France 7, 39 (1997). 
[8] P. Claudin, J. P. Bouchaud, M. E. Cates, and J. P. Wittmer, Phys. Rev. E 57, 4441 (1998). 
[9] V. M. Kenkre, J. E. Scott, A. Pease and A. J. Hurd, Phys. Rev. E 57, 5841 (1998) 
[10] S. Luding, J. Ceng, and R. P. Behringer, to be published. 

[11] J-P. Bouchaud, P. Claudin and D. Levine, M.Otto, preprint, cond-mat 0011213. 

[12] S. Luding, O. Straufi, in Granular Gases, T. Poschel, and S. Luding eds., Springer- Verlag, Berlin, 2000. 
[13] G. Reydellet and E. Clement, cond-mat/0012083. 

[14] J.J.Moreau, Colloque Physique et Mecanique des Milieux Granulaires (in french), Acte des Journees Scientifiques du LCPC, 

Sept 2000, pl99-204 (2000). 
[15] J. Rajchenbach, Nature 406, 708 (2000). 

[16] G. W. Baxter, R. P. Behringer, T. Fagert, and G. A. Johnson, Phys. Rev. Lett. 62, 2825 (1989). 



3 



20 



CI 

c 

ID 



X 

O. 



Fit Line Slope (a.u.m/N) 
□ (Dash Line) Pentagon^. 62 
• (Solid Line) Disk=5.45 




• Disk Exp. 
□ Pentagon Exp. 

- Disk Theory 

- Pentagon Theory 



20 



Depth (cm) 

FIG. 1. Hydrostatic pressure due to gravitational force alone versus depth. The expected slopes of the stress-height curves 
are calculated from the known packing fraction (7 = 0.91 for disks; 7 ~ 0.75 for pentagons). Inset shows the multi-particle G 2 
calibration by applying known loads to the upper surface of the layer. 
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FIG. 2. Images for pentagons showing a) the force pattern due to hydrostatic head, b) the combined response to gravity and 
a point source, c) the response to a point force-after subtracting image a from image b. d) similar to c, but with a different 
grain configuration. 
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FIG. 3. Mean response for 50 trials of a 50 g point force for (a) a uniform hexagonal packing of disks, (b) a bimodal packing 
disks with A = 0.99, and (c) for pentagons. The size of each image is about 18.0 x 13.5 cm 2 . 
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FIG. 4. Photoelastic response (force per length-stress in 2D) to a point force, vs. horizontal distance, x, at various depths, 
z, from the source, (a) for ordered disks, (b) for bimodal disks, and (c) for pentagons. Also shown are fits of the response to 
the convection-diffusion model. In all fits, c is dimensionless and D is in units of pixels (1 pixel = 0.45 mm). Inset part c: W 
vs. 2 for the response for pentagons. 
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